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STABILITY DERIVATIVE MEASUREMENTS WITH , 
MAGNETICALLY SUSPENDED CONE-CYLINDER MODELS 

D.Bharathan and S.S.Fisher + 

University of Virginia, Charlottesville, VA 22901 


ABSTRAC r 


In an initial feasibility study, the 
stability derivatives, c , C + C 

C 2 and C z> c for 5- and 7-caliber 

cone-cylinder models have been measured at 
M = 0.071 and Re = 1.3x10 4 by suspending 

each model electromaanetically in a small 
subsonic find tunnel , forcing it in 
periodic, combined pitchina/heavinq motion 
at freouencies near pitch resonance, and 
comparing its frequency response with flow 
to that without flow. Draq coefficients 
are measured as well. The apparatus and 
techniques employed are described, the 
analytical model used to extract the 
derivatives from the response data is 
outlined, typical response data are shown, 
comparisons are made with conventionally 
obtained, similar data from other 
facilities, and a general assessment of 
the technique is made. 


INTRODUCTION 


Preliminary stability evaluations for 
aerodynamic vehicles are usually carried 
out either with stinq-mounted models in 
wind tunnels or with models in free 
flight. In wind-tunnel tests, the 
influence of the sting can he significant 
and in free-lliqht tests observation time 
usually is severely limited and data 
retrieval is difficult. On the otherhand, 
magnetic suspension of wind-tunnel models 
eliminates the fluid-mechanics aspects of 
sting interference without incurring the 
complexities of free-fliqht testing. 

^The present investigation represents 
an initial application to stability 
testing of an electromagnetic suspension 
system which has been operated at the 
University of Virginia over the past ten 
years. It follows onlv a very few 
stability investigations of somewhat 
^limited scope conducted with magnetically 
suspended models elsewhere , -ft-4-) . -^These 
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4tests are aimed primarily at evaluating 
the practicality ot measuring stability 
derivatives with this apparatus. For the 
experiments, the nominal objective was to 
measure the pitchinq-momont derivative 
and, particularly, the pitch-dampinq 
derivative for each model. After an 
analytical model for the motion was 
developed and the data were obtained, it 
was found that side-force and side-moment¬ 
damping derivatives also could be 
estimated. 



The electromagnetic coil arrangement 
for the suspension system is shown jn 
Figure 1. A large pair of Helmholtz coils 
produces a uniform field to magnetize the 
model and a second pair of opposed coils 
creates a streamwise (vertical) gradient 
in the field at the configuration 
centroid, the nominal model-support 
location. This gradient induces a force on 
the model which opposes its weight and 
drag. Third and fourth sets of coils 
create lateral (horizontal) gradients in 
the field which induce side forces on the 
model. The Helmholtz coils carry up to 
200 A and produce a field of up to 120,000 
A/m at the model location. The vertical 
gradient coils carry up to 150 a and 

produce a field gradient as hiqh as 
390,000 A/m 2 at the model. The horizontal 
gradient coils carry up to 45 A and 

produce a gradient as hiqh as 9,500 A/ m 2 
at the model. All coils are water cooled. 







The wind-tunnel test section is 
inserted inside the group of coils. This 
tunnel is a small, fan-driven, open-return 
tunnel with hiqh contraction ratio. The 
test section is 11.4 cm b” 11.4 cm bv 30.0 
cm long, the air speed for these tests was 
24 m/s, and the flow turbulence level was 
0.4%. 


Model position is sensed optically 
using parallel light beams and silicon 
photocells. A feedback control svstem is 
then used to keep the model centered near 
its nominal support location. For this, 
feedback signals indicative of the model's 
position along each of three orthogonal 
axes, after passing through appropriate 
compensating electronic networks, drive 
power amplifiers supplying current to 
respective sets of force coils. 


For the stability tests, the model is 
forced periodically in combined pitching 
and heaving motion. Fere, model pitch 
orientation is monitored bv sensing its 
lateral position at separate stations 
along its length. To reduce the efrects 
of noise (due primarily to flow turbulence 
and secondarily to structural vibration 
and control-system noise), a 256-channel 
digital signal averager was used and the 
several position-sensor and coil-current 
waveforms were averaged simultaneously bv 
multiplexing the input to the averager. 


Other details concerning the 
apparatus as well as a more complete 
review of its operation and development 
are given in Ref. 5. 


MODELS 


The two models emploved are simple 
cone-cylinders. For each, the nose is a 
sharp cone 15.2 mm long with 14.5 degree 
semivertex angle. The cylindrical 
afterbody, 7.8 mm in diameter, is 22.9 nm 
long for the 5-caliber model and 38.3 mm 
long for the 7-caliber model. The cvlinder 
is hollow, with an inner diameter of 6.9 
mm, and is machined from a single Lucite 
rod. The magnetic part of each model is a 
hollow steel (A1SI-01) cvlinder 6.9 mm OD 
by 6.3 mm ID by 12.7 mm long which fits 
tightly inside the outer aerodynamic 
shell. For each model, the downstream end 
of the shell and the core are coincident 
and the model base is carefully covered 
with plastic tape. 


The models are made axisvmmetric 
because model roll is uncontrolled in this 
suspension svstem, and position control 
and sensing would be prohibitively 
difficult for a rolling, non-svmmetric 
model. The overall size of the model and 
core are limited by the wind tunnel's size 
(whicli is in turn limited by the space 
available between the coils) and the size 
of the region over which the magnetic 
field gradients are approximately 
constant. The relative sizing and mass of 
the aerodynamic shell and magnetic core, 
together with the imposed magnetic field, 
are chosen to place the model's pitch- 
resonance frecruency above the bandwidth of 
the position-control system, allow for a 
reasonably large aerodynamic drag, and 
keep the electromagnetic damping small. 
This selection leaves the control system's 
effectiveness essentially unaltered, while 
permitting adequate forced-motion 
amplitudes to be achieved rather readilv 
near pitch resonance. 


PROCEDURE 


For the tests, the model is driven 
magnetically with a sinusoidal lateral 
force over a series of freouencies 
covering pitch resonance. Due to the fact 
that the centers of mass ot the model 
shell and core are not coincident, the 
excitation force induces a combined 
pitching and heaving notion (predominantly 
pitching, since this ■ motion is lightly 
damped). To simplify interpretation of 
the overall motion, the pitch-angle 
amplitude is alwavs kept small, less than 
0.5 degree. Thus non-linear aerodynamic 
and electromagnetic effects on the model 
motion are minimized. This also allows a 
number of cross-coupled electromagnetic 
forces and moments acting on the model to 
be neglected. A theoretical model for the 
model's response is then used to extract 
estimates for unknown aerodynamic and 
electromagnetic coefficients. This model 
provides analytical expressions for the 
model's transfer functions in pitch and in 
heave. These expressions are algebraic 
functions of the aerodynamic and 
electromagnetic coefficients and the 
unknown coetficients are determined bv 
fitting the measured data to these 
expressions in a least sauare sense. 


Transfer-function amplitudes and 
phases are obtained upon computer 
processing the recorded averaged position¬ 
sensor waveforms and coil-current 
waveforms. In this, measured position¬ 
sensor and force-current calibrations are 
used to convert raw-signal data to 
positions, angles, velocities, and forces. 
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Data processing fields plots versus 
freguency of ( 1 ) transfer-function 
amplitudes, 6 o/Fo and w 0 /F 0 where e 0 is 
the pitch amplitude, w 0 is the non- 
dimensional heave velocity amplitude, and 
Fo is the non-dimensional perturbinq force 
amplitude, and ( 2 ) transfer-function phase 
angles and $ , the angles by which the 
pitch orientation and the heave velocity 
lead the perturbing force. * 

Transfer function variations are 
measured both with flow and without flow. 
Simple flow-oft (rather than vacuum) 
measurements are satisfactory here because 
still-air dampinq is negligible compared 
to electromagnetic damping. From the flow- 
off data, the unknown electromagnetic 
force and moment coefficients are 
determined and, once these are known, the 
aerodynamic force and moment coefricients 
can be determined from the flow-on data. 


Here, 

a = axial distance between the model's 
center of mass and center of mag¬ 
netization, 

C m = aerodynamic pitching 

a a q moment derivatives, 

*"z ,C z*' C z = aerodynamic side-force 
a a q derivatives, 

g = gravitational acceleration, 

I = vertical gradient-coil current, 
i B = model moment of inertia, 
k = magnetic pitch stiffness, 
n = core shape factor,(l-3N x )/2(1-N x ) 
N x = core axial demagnetizing factor, 
s = Laplace transform variable, 

Y x = vertical coil force constant, 
e m = m f 9 net ic damping moment coeffi- 
e = magnetic damping force coeffi- 
m = model mass. 


THEORETICAL TRANSFER FUNCTIONS 


In Ref. 5, the following analytical 
expressions for the model transfer fun¬ 
ctions are derived: 


where 


w Pj I s 2 + 2 P 2 P 3 S + P§ ) 

F s [ s 2 + 2PsPss + Pi ] 

0 P4 

F [ s 2 + 2P s P 6 s + Pi] 


(1) 

(2) 


Pi = 1/2)1 

Pi = I k + 2uga - C m - aC., ]/i Q 


2P 2 P: 


2P 5 P 6 


(C +C )-(C in .+C m )]/i B 
d z q a q 


[k+Y x I XJ (l+n)a - C m ]/i B 
[e m +a 2 c z -(C m> +C m >]/i B 


In these expressions, P 3 and P 6 rep¬ 
resent natural frequencies, P 2 and P 5 
represent the associated damping factors, 
and Pi and P, are proportionality cons- 


* In this work, velocity is normalized 
with the airspeed, V, and electromagnetic 
lateral force is normalized with the 
flow dynamic pressure times the model 
base area. For the aerodynamic parameters, 
the model base diameter is used as the 
normalizing length. 


In developing these expressions, as 
justified in Ref. 5, the effect of the 
heaving, motion on the pitching motion is 
neglected. i.e. the pitching motion is 
treated as that for a simple linear damped 
oscillator. Also, magnetic dampinq in 
pitch and heave are assumed linearly 
proportional to the pitching anqular 
velocity and the heaving velocity, 
respectively. For the small motion 
amplitudes here, these latter assumption 
are satisfactorily accurate. 


TRANSFER FUNCTION DATA 


Measured flow-off variations for the 
pitch transfer function for the 5-caliber 
model are shown in Fiqure 2. The data in 
this figure and following figures are 
taken from four sets of independent runs, 
two each for excitation alonq a separate 
lateral axis. Since within the scatter in 
the data the four sets identical, all data 
have been combined. For presentation, the 
phase data is plotted as cos <)> rather than 
<l> itself. The plotted curve through the 
amplitude data is a least-sguares fitted 
version of Eq.(2). A standard computer 
code ( 6 ) was used to obtain this fit. 
Using the values of Pj (1*1,6) determined 
for this match, the curve through the 
phase data is obtained, i.e. the phase 
curve is not a least-sguares fit to the 
phase data , but is that corresponding to 
the P • values determined from the 
amplitude match. Of course, the plotted 
curves in this figure are merely those for 
resonance ot a simple damped oscillator. 
The amplitude exhibits a sharp peak at the 
resonant frequency and across this peak 
the phase chanqes rapidly by 180 deqrees. 
The excellent agreement between the 
plotted curves and the data is convincing 
evidence that the pitching motion is 
properly modeled. 



































STABILITY DERIVATIVES 


and 2. Values of the axial force (neqative 
drag) coefficient, C x t also measured in 
these experiments, are listed as well. 
Analytical estimates and comparable 
measurements obtained elsewhere are also 
included. The analvtical estimates are 
taken from the t'SAF DATCOM (Data 
Compendium) (7) and, except for C x are 
based on an assumption of inviscid flow 
over a slender body. Mo comparable data 
for C +C_ for either model could be 


5-caliber model. Also, : 
data jure for exact ma 
model shape. These 
reasonably close, hovev 


uncertainty, is 15% below its analytical 
estimate and 10% below the mean of the 
other measurements. The C measurement, 
with 2% uncertainty, lies between the two 
analytical estimates, one assuming fully 
turbulent and the other fully laminar 
flow. This suggests a partiallv turbulent 
boundary layer over this lonqer model. The 
C + C measurement, with 25% 


with essentia)ly 100% uncertainty, is 
of its analytical estimate. As alr< 
mentioned, no comparable data for c 


value, with 60% ui 
uncertainty), is 75% of tl 
measurements and 20% belo\ 


uncertainty, is 3% above 
measurements and 10% abc 
estimate. The c x Vi 
uncertainty, compares 
laminar analvtical 


The present 
the feasibility 
stability deriv; 
They also il] 
complexities in\ 


r analvtical 
ting a fully 1, 
this model. 
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5-Caliber Model 












errors in the model ari 
than those correspondii 
the response data. 
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